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ABSTRACT: 


We  consider  convergence  of  an  approximation  method  for  the  recovery 

of  a  rotationally  symmetric  potential  ^  from  the  sequence  of 

eigenvalues.  In  order  to  permit  the  consideration  of  'rough'  potentials 

'P  (having  essentially  H  *(0,1)  regularity),  we  first  indicate  the 

appropriate  interpretation  of  -A-h^  (with  boundary  conditions)  as  a 

2 

self-adjoint,  densely  defined  operator  on  3f  :=  L  (fi)  and  then  show  a 
suitable  continuous  dependence  on  'p  for  the  relevant  eigenvalues.  The 
approach  to  the  inverse  problem  is  by  the  method  of  'generalized 
interpolation'  and,  assuming  uniqueness,  it  is  shown  that  one  has 
convergence  to  the  correct  potential  'p  (strongly,  for  an  appropriate 
norm)  for  a  sequence  of  computationally  implementable  approximations 


2 


1  INTRODUCTION 

The  present  paper  is  intended  as  an  extension  of  the  considerations 
of  [5]  to  higher  dimensional  contexts.  Our  concern  will  be  with 
operators  of  the  form 

(1.1.)  L  =  L^:u-^-V*aVu  +  *u 

in  a  context  of  rotational  symmetry  in  ,  i.e.  assuming  that 
(1.2.)  (i)  a (*)•  'M*)  depend  only  on  r  :  =  |x|, 

(ii)  the  domain  12  is  the  unit  ball  of  IR^  with  d  >  2, 

(iii)  the  boundary  conditions  are  radial  -  of  the  form^ 

au^  =  tu  on  3  !)  (i.e.  at  r  =  1)  . 

We  assume  a(*)  is  known  and  bounded  with  a  uniform  ellipticity 
condition: 

(1.3.)  A  >  a(r)  >  a  >  0  for  0  <  r  £  1 

(e.g.  a  =  1  giving  Lq  =  -  A) . 

Our  concern  is  with  the  inverse  eigenvalue  problem.: 

(EVP)  Suppose  a(*)  is  known  and  it  is  known  that  'p  €  (some 
suitable  set).  If  we  are  given  eigenvalues  of  the 
self-adjoint  operator  associated  with  (1.1)  and 

(1.2)(iii),  how  can  we  (computationally)  recover  the  potential 
+  ? 


^We  could  equally  well  consider  Dirichlet  conditions  (u  =  0  at  r  =  1) 
which  would  require  minor  modification  of  our  presentation  -  e.g.  we 

would  set  V  :  =  Hq  (12)  rather  than  H^(12)  as  here,  etc. 
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We  do  not  consider  here  the  deep  question  of  uniqueness:  within  which 
sets  P  is  P  uniquely  determined  by  the  given  eigenvalue 
information.  Rather,  this  is  taken  as  an  a  priori  assumption  on  the 
suitability  of  P*  for  (EVP). 

On  the  other  hand,  as  in  [5],  we  are  very  much  concerned  with 

another  aspect  of  the  suitability  of  P  :  For  how  ’rough’  a  potential 

\ p  can  we  construct  a  workable  interpretation  of  (-  V  *  a  V  +  p)  as  a 

2 

densely  defined,  self-adjoint  operator  on  L  (Q)  with  compact 

resolvent  so  that  discussion  of  the  "eigenvalues  of  ...  A^”  makes  sense? 

The  approach,  as  in  [5],  is  closely  related  to  that  of  Chapter  3  of  [6] 

with  modification  to  fit  the  setting  under  consideration.  In  [5]  it  was 

shown,  for  the  one-dimensional  case,  that  A^  is  suitably  defined  for 
^  1 

p  €  S6  with  9s:  =  H  (-1,1)  and  that  the  eigenvalues  then  depend 

2 

continuously  on  p  .  A  principal  concern  here  will  be  to  obtain 

**  1 

comparable  results  for  p  €  9  with  9s  much  like  H  (0,1)  -  viewing 
p  =  ^(r)  as  given  for  r  €  (0,1),  rather  than  on  0  -  but  now  with  9 
defined  through  a  weighted  -  norm,  controlling  the  behaviour  near 
r  =  0  .  We  are  able  to  get  results  quite  comparable  to  the 
one-dimensional  case  treated  in  [5]  precisely  because  the  radial 
symmetry  permits  a  treatment  through  separation  of  variables  which 
reduces  this  to  one-dimensional  considerations. 

Once  we  have  developed  the  setting  in  which  (EVP)  is  a  meaningful 
problem,  our  concern  is  to  demonstrate  convergence  for  an  approximation 
method  of  ’generalized  interpolation’  type  (see,  e.g.,  [4], [6]).  We 
assume  in  (EVP)  that  we  are  given  the  sequence  (A^.Ag..--)  of  the 
eigenvalues  of  A^j  corresponding  to  a  unique  potential  p  €  P*  and 
consider  the  approximation  procedure: 


^There  appears  to  be  a  slight  gap  in  the  argument  in  [5]  and  it  seems 
necessary  to  take  the  strong  topology  for  P  rather  than  the 

sequential  weak  topology  as  asserted  there.  See  Theorem  11,  below. 
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(ii)  For  a  €  L°°  and  ^  in  a  certain  Lq(Q)  (see 
below)  we  have  a  'weak  interpretation'  of  (L^+X)  as 

a  continuous  invertible  operator  ■  V  — *11  with 
II  =  H^n)  .  We  can  then  take  to  be  the 

pre-image  of  %  C  f  for  this  operator. 

(iii)  The  radial  nature  of  induces  (e.g. ,  for  smooth 

a.vp  )  a  canonical  decomposition  of  3f  by  separation 
of  variables  into  subspaces  of  the  form  H  =  3CWI 

H"  r1 

3 

where  each  y  is  finite-dimensional  and  9L  is  a 
M 
o 

weighted  L  -  space  of  functions  on  (0,1). 
Associated  with  -V*aV»  and  (1.2)(iii)  is  an 
ordinary  differential  operator  M  and,  following 
[5],  we  can  interpret  (M^-h//)  as  a  self-adjoint 

operator  on  2C  for  each  relevant  p  when  ^  is  in 

^  1 

where,  now,  &  is  a  weighted  H  -  space  on 

(0,1).  These  interpretations 

A  :  K  D  2J  ,  — *  If  can  be  combined  to  obtain 
-M.'#'  M  M.'/'  M 

SD,  =  2)  (***)  and  the  interpretation  of  (2.1). 

'P  Y 

We  will  ultimately  use  the  interpretation  (iii)  but,  of  course,  wish  to 

know  that  the  interpretations  are  consistent  with  each  other. 

We  begin  with  an  abstract  construction,  following  [5].  For  the 

2 

present,  we  take  3t  to  be  L  (if)  where  if  is  any  set  with  (positive) 

M 

bounded  measure  n  and  assume  *3/  is  also  a  Hilbert  space  of  functions 
on  if  with  the  pivoting 

y  S£  <*  y* 


3 

The  elements  of  II  are  just  the  classical  'angular'  functions 

well-known  from  analysis  of  the  Laplacian  for  a  ball  -  i.e., 

(sin  n0,cos  n0}  for  d  =  2.  spherical  harmonics  for  d  =  3,  etc.  The 
subspaces  {*31  }  do  not  depend  on  a , 41  or  the  boundary  conditions  and 


give  an  orthogonal  direct  sum  decomposition  of 


L2(an)  . 


with  dense  embeddings.  We  assume  given  a  linear  continuous  map 
M:  <3/  — *  y  for  which  one  has  a  monotonicity  estimate  of  the  form: 

(2.3)  <  Mx,  x  >  ^  a  |x | J  -  0  |x|^  x  €  V  C  2C 
with  a  >  0  .  We  will  also  assume  symmetry: 

(2.4)  <  Mx.y  >  =  <  x,My  >  for  x.y  €  <2/  . 

(The  <2/  -  *2/*  dualities  of  (2.3),  (2.4)  are,  of  course,  given  by  the 
pivoting  through  the  inner  product  of  3 C  ). 

We  next  take  9  to  be  any  space  containing  products  xy  for 

4 

x.y  €  <2/  with  a  norm  such  that 

(2-5)  |xy|^  <  Clxl^lyl^  for  x.y  €  <2/  . 

LEMMA  1  Let  ^  be  a  function  on  9  which  is  in  9  in  the  sense  of 

x 

the  9-9  duality  induced  by  the  9C  inner  product.  Then  the 
multiplication  operator 

x  -*  \px'-  y  -+ 

is  well  defined  and  continuous  with 

(2.6)  N*l„*  1  Cl^lxl^  so  il^ll  S  Cl*^*  - 

Proof:  From  (2.5)  we  have 

K'//x,y>  I  =  I  <'p,xy>  I  £  l'Pl<p*  Ixy  1^  ^  (C  |x|/ 1^*  |x  |^)  ly  1^ 

for  arbitrary  x.y  €  <3/  .  By  the  definition  of  the  <2/  norm  as 
sup{  |<v//x,y>  | :  |y  1^  i  1}  ,  this  gives  (2.6).  □ 

^E.g.  We  could  take  9  to  be  the  completion  of  sp{xy:x,y  €  <2/}  with 
respect  to  the  norm 
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Note  that  for  such  functions  \p, <p  we  have  <f  >  ^  precisely  when 

9 

y.y>  =  <'P~'P . y  >  >  o  for  y  €  <3/  and  we  take  this  as  inducing  the 

order  for  #  .  We  also  wish  to  consider  <f  €  ^  (so  that  <g:*3/  -»  <3/  is 

defined)  such  that: 

(2.7)  For  each  e  >  0  one  has  C£  such  that 

|<<px,x>|  £  e|x|^  +  C£|x|^  for  x  €  <3/  C  91  . 

LEMMA  2  Let  3C , <3/ . ^  be  as  above.  Let  M:<3/  -»'?/*,  as  aboue,  satisfy 

(2.3),  (2. 4)  and  let  \p  €  $>*  with  2  <p  for  <p  satisfying  (2.7). 

Then  (M  +  \p)  induces  a  densely  defined,  self-adjoint  operator 
M^:91  D  2)^  -»  91  .  If  the  embedding  <5/  «*  21  induced  by  the  pivoting  is 

compact,  then  has  compact  resolvent . 

Proof:  For  any  real  A  we  have  (A+vJ>+A)  :<3/  -*  “3/  continuous  with 

2 

<(M-H/rf-A)x,x>  =  <Mx,x>  +  <A+<#>)x,x>  +  <^p->p,x  > 

>  [glxlj  ~  Plxlj]  +  Mxlg  +  <wx.x>  , 

using  (2.3)  and  noting  'P  >  f  ■  Using  an  inequality  |x|<£  £  Co^x^<2/ 

(2.7)  with  e:  =  a/2c  ,  we  obtain  the  fundamental  estimate 

v  *  ~  o 

(2.8)  <(M+^A)x,x>  l  (a/2)  !x|^  +  (A-0-^C. )  |x  . 

—  ^ 

Considering  A  >  /3  +  cqC&  =■  •  this  makes  (M-hJh-A)  :  V  -*  9/  strictly 

monotone,  hence  invertible.  We  set 

(2.9)  ^:=  a((M+*+A)_1|a) 

=  (x  €  V:  (M+\^+A)x  =:  z  €  1}  Cl  , 

M^x:=  z-Ax  for  x  €  2)^  with  z:=  (M+v|/+A)x  €  91 
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_  1  ^ 

The  continuity  of  (A+v/h-A)  :  X  *»  2/ 


y  ensures  that  A,  is  a  closed 

~\U 


operator.  Clearly,  this  definition  is  independent  of  the  particular 
choice  of  (large  enough)  A  .  Note  that  A^  depends  on  only 
through  the  lower  bound  <p  but  (2.9)  does  not  depend  on  the  particular 
choice  of  <p.  Fixing  <p  and  taking  A  >  A  ,  (2.8)  gives 


(2.10) 


Ixl^  ^  la  for  x  €  2^ 


uniformly  on  (^  €  &  '■'p  >  <p}  with  M  =  2/aC  where  C:  =  [norm  of 
embedding:  2/  -»  3 C]  .  Note  that  if  2/  -*  X  is  compact  then  (2.10)  makes 

~*1  _  H 

{(M^+A)  ;A>A^,'//>ip,\^e9i}  collectively  compact  for  any  <p 

satisfying  (2.7). 

To  see  that  2)  is  dense  in  X  ,  i.e.  that  (M+^h-A)  * :  X  -»  X  has 
dense  range,  we  proceed  by  contradiction.  Were  2 3  not  dense  there 
would  exist  x  €  X  orthogonal  to  3^  with  x  ^  0  .  We  can  find  z  €  2/ 
with  (M-h/h-A)z  =  x  so,  using  (2.4), 

<z,  z>  =  <(M+v/rt-A)z,  (M-H/rt-A)_1z> 

=  <x,  (M+^+A)  ^z>  =  0  since  (M+^h-A)  *z  €  2)^ 
for  any  z  €  X  .  Hence  z=0  so  x=0-  a  contradiction. 


Finally,  the  assumed  symmetry  of  M  makes  formally 

self-adjoint  but  we  must  verify  that  the  domain  of  (M^)  is  precisely 
2!^  ,  i.e.  that  X-continuity  (on  the  dense  set  2)  )  of  the  functional: 
y  -*  <M^y,x>  implies  x  €  2)  (noting  that  the  inverse  implication  is 
clear).  Note  that  this  continuity  implies  existence  of  z  €  X  such 
that  <M^y,x>  =  <y,z>  for  each  y  €  2 1^  ;  set 

x  :  =  (M-h/h-A)_1[z+Ax]  €  ffl  . 


We  then  have,  by  the  symmetry, 
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<(M+'/'+A)y,x>  =  <My,x>  +  <('//+X)y,x> 
=  <y,Mx>  +  <y,('|H-X)x> 
=  <y,  (M+^X)x> 


=  <y,z+Xx>  =  <M^y,x>  +  X<y,x> 
=  <(M+^X)y,x>  . 


Since  (M+v/n-X)y  ranges  over  1  as  y  ranges  over  2)^  ,  this  shows 

x  =  x  so  x  €  2),  .  0 

'P 

cl  1 

For  the  interpretation  (2.2)  (ii)  of  ,  we  take  °y  =  If  :=  H  (Q) 

2 

and  3t  =  It  ■  =  L  (Q)  .  The  standard  weak  formulation  of  A  :=  -  V-aV* 
with  (1.2)  (iii)  is  given  by 


(2.11) 


<Au,v>  =  I  (-V*aVu)v 

=  aVu*Vv  -  -r  v 

Jn  Jan 


Since  [1]  the  Dirichlet  trace  is  compact  :  <?/  -*  Lr(90)  ,  we  have 


(2.12) 


ij  ■ 

J3Q 


uv|  i  elul^lvl^  +  ceiul3flvlgf  for  u>v  €  V 


by  a  standard  functional  analysis  result.  Hence,  taking 


5  2 

As  noted  earlier,  we  have  H  ■  =  Hq(Q)  in  the  case  of  Dirichlet 

boundary  conditions. 

0 

Given  Banach  spaces  %,1I ,f  and  linear  maps  A  :  °U  -*  II  and  B  :  W  -»  W 
with  A  compact  and  jV(B)  C  .Af(A) ,  one  has  an  inequality: 

|Au|^  <;  elul^  +  C£|Bu|^  for  u  € 

for  arbitrary  e  >  0  and  some  C£  . 


10 


(2.13)  luL  :=  [f  lVu|2  +  |u|2]1/2  . 


we  obtain  (2.3)  for  M  =  A  with  a  arbitrarily  close  to  the  a  in 

7 

(1.3)  and  a  correspondingly  determined  /3  .  We  now  take 
q  :=  2d/(d+2)  >  1  so  that  (with  1/q  +  1/p  =  1)  standard  results  [1] 

give  continuous  embedding  :  H  L  P(Q)  whence 

lx2 I  -  =  lx  I2  i  Clxl2  . 

P  2?  f 

This  means  that  we  can  take  $  =  LP(H)  and  \p  €  &  =  Lq(0)  .  For 
q  >  q  (or  q  =  q  when  d  =  2)  the  embedding  :  H  -*  L2p(Q)  (with  1/q  + 
1/p  =  1)  is  compact  so  we  have  (2.7),  for  <f  €  Lq(Q)  .  We  have  thus 
shown: 


(2.14)  For  any  ^  €  Lq(Q)  ,  i.e.,  Jr  1 |^(r) lqdr  <  ®  ,  bounded 

below  by  <p  €  Lq(n).  the  construction  (2.2)  (ii)  via 

(2.9)  defines  ,  corresponding  to  in  (1.1),  as  a 

2 

densely  defined,  self-adjoint  operator  on  L  (Q)  with 
compact  resolvent.  □ 


When  a,v//  are  smooth  it  is  a  standard  regularity  result  that  u  € 
for  u  €  L2(Q)  so  the  definitions  (2.2)  (i),  (ii)  are  then 
equivalent . 

g 

We  proceed  now  to  develop  (2.2)  (iii).  A  formal  calculation 
imposing  the  ansatz 


H2(Q) 


(2.15)  u(x)  =  R(r)U(w)  for  x  =  ru  €  Q 

with  r  €  (0,1)  and  co  €  ^  =  5ft  ,  gives 


7This  is  for  d  >  2  .  For  d  =  1  we  could  take  9  =  as  in  [5]  while 
for  d  =  2  we  take  any  <p  >  1  and  continue. 

Q 

This  is  valid  pointwise  for  a,R,U  smooth. 


11 


(2.16)  -  V*aVu  =  (MqR)U  +  -^R(SU) 

r 

where  is  the  ordinary  differential  operator  given  formally  on  (0,1) 

by 

(2.17)  Mq:  f  --  (af1)*  -  af * 

and  S  is  a  second  order  elliptic  operator,  acting  as  a  densely 

2  d-1 

defined,  self-adjoint,  semi-definite  operator  on  :=  L  (S  )  with 
compact  resolvent. 

A  significant  observation  is  that  the  spherical  operator  S  does 
not  depend  on  A  ,  i.e.  on  the  particular  choice  of  a(*),  nr  in  (1.1) 
and  (1.2)  (iii).  We  write  {/i^  :  j  =0,1,  ...}  for  the  distinct 
eigenvalues  of  Mq  so 

0  =  <  ...  -*  00 

and,  for  each  p  =  p .  ,  we  let 

J 

(2.18)  H  :=  (U:  SU  =  p  U}  C  11  :=  L2(Sd_1) 

be  the  corresponding  eigenspace.  Note  that  each  It  is 
finite-dimensional  and  that  the  elements  of  11 (eigenfunctions  of  S  ) 
are  just  the  classical  'angular'  functions.  The  subspaces  (11  }  are 
othogonal,  giving  a  direct  sum  decomposition: 


(2.19) 


^  =  {^:  M  =  0,p1 - }  . 


Now  let  J  be  a  weighted  L^(0,1),  i 


(2.20) 


r  1 


<x,y>^  := 


0 


d-1  , 

r  x(r 


and  corresponding  norm.  For  p  =  0,  p^ , . . 


.e.,  with  the  inner  product 

)y(r)dr 

we  set 
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(2.21) 


:  =  X&V  (tensor  product) 


:=  sp  (u  =  R(r)U(w):  R  €  2  .  U  €  <U  } 

r1 

r  j(m)  _  I  2 

=  il  R.(r)  U.(«):  R.  £1  C3M=  IT(n) 

Lj-i  j  j  j  J 

where  (U.  =  U  j=l . J(p)}  is  3X1  othonormal  basis  for  V  . 

\  j  [j  j  o  yr-jj  p 

Corresponding  to  (2.19)  we  then  have  an  orthogonal  direct  sum 


decomposition: 


(2.22) 


#  = 


M:  M  =  O.Mj - }  ■ 


Note  that  the  norm  |  ♦  corresponding  to  (2.20)  gives  |RU|^  =  |R|^|U|^ 
for  (2.15)  and  if,  corresponding  to  (2.-21),  (2.22),  we  consider  u,v  € 
expanded  as 

„  J(m)  _  v  . 

U  =  ^  Z..l  'V  =  ^ 


wi  th  R  . ,  ft  .  €  3C  ,  then 
M.J  M.J 


J(M) 


<u.v  <Vr 

in  view  of  the  orthonormality  of  (U  :  j  =  1 . J(p);p  =  0,/i.  ,...}  . 

M*  J 

Note  that  for  u  as  in  (2.15)  with  U  £  1  (e.g.  U  =  U  .)  we  have 

fJ.  J 

=  [(Mo+Ma/r2-Hl/)R]U 
2 

and  we  are  led  to  analyze  :=  Mg+pa/r  (p  =  p^.p^,...). 

Integrating  by  parts  and  using  the  boundary  conditions  (1.2)  (iii) 
we  obtain  the  weak  formulation  of  M  : 

~p 


<M 


f.g>».  =  f  rd_1a(r)[f 'g'+^-  fg]dr  -Trf(l)g(l)  . 
1  a  J0  T 


(2.23) 
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To  proceed  it  is  necessary  to  distinguish  the  two  cases:  pi  =  Mq  =  0 

and  pi  =  pt ^.pi^,...  >0  .  In  each  case  we  take  <3/  to  be  a  weighted 

space,  but  use  slightly  different  norms.  We  take  to  be  the 

9 

Hilbert  spaces  of  functions  f  on  (0.1)  induced  respectively,  by  the 
norms 


(2.24) 


llfll0  :=  ( 


d-lr  ...  ,2  ...2-, ,  ,1/2 
r  [If  I  +|f|  ]dr ) 


nr..  •  /•  f  1r  d-1 ...  .2,  d-3 . r  ,2-,  ,  ,1/2 

Ilf  il  :=  (  [r  If  !  +pi..r  If  I  ]dr) 

J0 

_2 

Observe  that  IMI+  dominates  IMIq  since  r  >1  and  pi^>  0  so 

<3/+  C  <8/  with,  clearly,  a  dense  embedding.  We  complete  the  weak 

formulation  of 


M^:  «->«**  (<3/  :=  for  pi  =  0;  V  =  <3/+  for  pi  >  0) 

by  specifying  that  (2.23)  is  to  hold  for  f,g  £  H  ,  as  appropriate. 

From  standard  embedding  results  [1],  one  easily  sees  that  <3/^  (a 

1/2  -  — 

fortiori  <3/+)  embeds  in  C  [r.l]  for  any  r  >  0  so  one  has  an 
estimate 


If (1)  I  <  C1IIf l!0  for  f  €  <3/Q  . 

Also,  *3 /  embeds  compactly  in  C[r,l]  for  r  >  0  from  which  it 
follows,  as  earlier  for  (2.12),  that 


9 

We  take  the  closure  (with  respect  to  IMIq  or  IMI+)  of  the  set  of 

smooth  functions  compactly  supported  in  (0,1]  and  satisfying  the 
boundary  conditions.  Hence,  the  specification  of  %  includes  the 
requirement  that  f(l)  =  0  in  the  case  of  Dirichlet  conditions.  Note 

that  for  pt  =  0  we  have  11^  =  {constants  on  *}  so  Kq  is  just  the 

2  ~ 

subspace  of  radial  functions  in  L  (Q)  while  (2.24)  makes  i3/q=‘3/q$<3 l q 

the  subspace  of  "f  '■=  H^(Q)  (or  Hq(Q))  consisting  of  radial  functions 
with  a  norm  isometry. 
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(2.25)  If  (1)  |2  <  e  Ilf  II2  +  C£  If  |2  for  f  €  VQ 

for  any  e  >  0  .  We  will  need  more  precise  information  about  the 
behaviour  of  f  €  ^  as  r  -»  0+  .  For  0  <  r  <  1  we  have 

If (1)  -  f (r)  |  =  |[  1s”(d"1)/2[s(d_1)/2f,(s)]ds| 

-  ,  P  *  1-d  ,  x  1/2  ..... 

i  (  s  ds)  IlfllQ 
Jr 

whence ,  as  |f(r)|  <,  |f(l)|  +  |f(l)  -  f ( r )  I  .  we  have 

1/2 

[C  +  in  r]IIf  ll0  d  =  2 

(2.26)  If (r)  I  £  • 

C H  r1-d/2llfll0  d  =  3,4 _ 

with  C^.C^  depending  only  on  d  .  If  we  define  as  the  space  of 

functions  continuous  on  (0,1]  for  which  the  norm 

(2.27)  *f*[>]  :=  suP{rU/2If(r)  I  :  0  <  r  i  1} 

is  finite,  then  (2.26)  shows  that  *3/q  embeds  (continuously,  by  the 
Closed  Graph  Theorem)  in  2^  for  u£i>:=d-2(u>u  =  0  for 
d  =  2)  . 

LEMMA  3  '5/q  embeds  compactly  in  for  v  >  v  :  =  d  -  2  (d  >  2)  . 

Proof:  Suppose  d  >  2  and  {f^}  is  bounded  in  .  We  can  extract  a 
subsequence  (again  denoted  by  {f^})  converging  “S/q  -  weakly,  say  to 
f  ,  and  we  will  show  f ^  -»  f  in  ^  .  Note  that  f^  ->  f  uniformly  on 
[r,l]  for  each  r  >  0  since  the  embedding  :  “3/q  -»C[r,l]  is  compact 
(as  II * IIq  dominates  the  H^(r,l)  -  norm);  cf.  e.g.  [1].  Since  v  >  v 
we  have,  by  (2.26). 
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(2.28)  r”/2|fk(r)  -  f (r)  |  1  C*  rU’llfk  -  fllQ  £  CrU' 

with  2d  1  =  v  -  v  >  0  and  with  C  f ixed  f or  the  sequence .  Given  any 
e  >  0  we  can  choose  r  so  the  right  side  of  (2.28)  is  less  than  e  on 
(0,r)  .  Then,  noting  the  uniform  convergence  fk  -»  f  on  [r,l]  ,  we 
can  choose  K  =  K(e)  large  enough  that  the  left  side  of  (2.28)  is  less 
than  e  on  [r,  1]  for  each  k  >  K  ,  giving  |fk  -  £  e  . 

<v 

For  d  =  2  we  maytake  any  u  e  (O.d)  and  note  that  (2.26)  gives 

IV 

If  (r)  I  <  r“u/2llfll0 

(C^  now  depending  on  the  choice  of  u),  giving  (2.28)  with 

IV 

2d'  =  d  -  d  .  The  proof  concludes  as  before.  0 

Next,  we  consider  the  weighted  (0,1)  spaces  induced  by 

the  norms 

P  1 

(2.29)  llfll  :=  (  rU[|f' |2  +  If  |2]dr)1/2 

L  J  Jq 

for  d  >  0;  observe  that  =  #  .  We  will  set  :=  <!/_  with 

v 

v  •=  2d-3  for  d  >  2  (any  d  >  1  for  d  =  2)  . 

LEMMA  4  Let  <p  €  (#  )  for  some  v  >  2d-3  (d  >  2). 

Then  (2.7)  holds  with  y  =  3/q  (a  fortiori  with  ‘J/  =  <3/+). 

Proof:  Set  v  ■=  v  -  (d— 1 )  so  v  >  2d-3  ,  as  assumed,  gives  d  >  d-2 
whence,  by  Lemma  3,  the  embedding  :  *3/q  Z~  is  compact.  Again  this 
gives  an  estimate 

IV 

(2.30)  rV  If (r)  |2  £  e  llfll2  +  C  |f|2  for  f  e  C  3E 

for  arbitrary  e  >  0  and  with  C  depending  on  e,  v,  d  . 

(f2)'  =  2ff '  so 


We  have 
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JrU  I  (f2) '  |2dr  £  4  sup{rU If (r) |2}  Jrd  1|f'(r)|2dr 
i  4  (anfn2+eif  i2)  if  i2 

rU|f2(r)|2dr  i  (£llf  ll2+ft  |  f  |2)  I  f  |2  . 

Ilf2ll2u;i  *  4(£llfll2+ft|f  |2)llfll2  . 

Setting  4e  :=  e  (also  fixing  ft)  ,  we  then  have 

(2.31)  llf2ll[u]  £  e  llfll2  +  C|f  |2 

for  any  C  large  enough  that  2eC  £  e  +  4ft  and  also  C2  >  ft  . 

For  <p  €  (“2/^)  let  ^  :=  norm  of  <p]  <  00  and  note  that 

(2.31)  gives 

|<<px,x>  I  =  l<*,x2>|  ^  Mllx2llj-u-j  <(eM)|x|2+  (CM)lx|2 

for  x  €  <2/  =  (or  <3/ +C  <3/^)  and  any  e  M  >  0  and  correspondingly 
determined  CM  .  This  is  just  (2.7).  0 

LEMMA  5  With  9  :=  ,  as  above,  we  have  (2.5)  rnith  <2/  =  *2/  (a 

v 

Fortiori,  with  <2/  =  *3/  ) . 

Proof:  Note  that  v  '■=  v  -  (d-1)  >0  for  any  d  >  2  .  We  no  longer 
have  compactness  but  (2.26)  gives 

rU|f(r)|2  <  Cjf  |2  for  f  €  <2/ 

corresponding  to  (2.30)  so  we  obtain,  as  for  (2.31),  the  estimate 
|x2|^  :=llx2ll|-— -j  £  C*|x|2  for  x  €  <2/  , 
absorbing  Ix terms  in  Ixl^  •  From  the  identity 

xy  =  [(cx+y/c)2  -  (cx-y/c)2]/4 


17 


we  then  obtain 

Ixyl  1  (C*/4)[ Icx+y/c 1^  +  | cx-y/c I 

£  (C*/2)[c  |x  |^+  |y  l^c]2 

_  ^  2 

which  is  just  (2.5)  with  C  :=  2C  on  setting  c  :=  lyl^lxl^  .  □ 

With  these  lemmas  in  hand  we  are  ready  to  proceed  to  the 
construction  (2.2)  (iii). 

THEOREM  6  Let  a(*)  satisfy  (1.3)  with  the  boundary  cond.it ions  (1.2) 

/V 

(iii).  With  ty  :=  %_  as  above,  assume  'p  €  S6  with  'p  >  <p  for  some 

v 

as  in  Lemma  4.  Then  L^  =  A  +  ^  induces  a  closed,  densely  defined, 

self-adjoint  operator  on  3f,  as  in  (2.1). 

Proof:  We  begin  by  considering  the  ordinary  differential  operator 
(M^-h£)  ,  first  for  p  >  0  so  we  are  taking  %  =  V+  and  I  •  1^  =  IMI+  . 
From  (2.23),  using  (1.3)  and  (2.25),  we  have 

<M  x,x>  >  a  f  [r^  ^  |x' +  (J.  ^|x|^]dr  -  t+|x(1)I^ 

Jq 

>  alxl^  +  (p-M1)alx|^  -  T+[e|x|J  +  Cjxl^] 

where  n  :=  max  {-r.O}  .  Choose  e  <  a/~r  (thus  also  fixing  C  )  and 
one  obtains  (for  p  >  1)  the  monotonicity  estimate^ 

(2.32)  <M^  x.  x>  ^  a  |x|^  +  (op  -  P0)|x|^  for  x  €  *  =  «J+ 

where  a  :=  a  -  and  :=  +  -r+C£  .  i.e.  (2.3)  holds.  The 

symmetry  condition  (2.4)  is  clear  from  (2.23)  and  Lemmas  3  and  4  ensure 

^We  only  need  (2.3)  immediately  but  emphasise  that  a  and  /?q  are 
independent  of  p  for  n  =  . . .  . 
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the  hypotheses  on  >//,  <p  for  applicability  of  Theorem  2.  Thus  we  know 
that  for  each  n  =  i±^,  ■  .  ■  there  is  a  well-defined  self-adjoint 

operator 


(2.33) 


m  .  ms  .  a 


which  maps  :  x  -»  (Mx  +  \//x)€  a  whenever  x  €  2)  C  ?  C  I  .  For 

jj.  %\jj  + 

fj.  =  Pq  =  0  ,  taking  2/  =  2/q  and  1*1^  =  H*Hq  •  one  similarly  obtains 
(2.3)  with  a  and  £5  =  J3q  exactly  as  in  (2.32).  Thus,  we  have  (2.33) 
for  euery  p  =  Mq.M1 . 

At  this  point  we  recall  (2.21),  (2.22)  and  note  that  each  u  € 
has  the  orthogonal  expansion 

(2.34)  (*  =  ™) 


where  the  outer  sum  is  over  p  =  lu,  u. , . . .  and  each  R  .  is  in  a 

0  HI  M.J 

The  orthonormality  of  (U  .}  gives  the  norm  identity 

M » J 


(2.35) 


lu|2  =  X  IR  .1 

H  H=1  M.J 


with  a  corresponding  formula  for  the  inner  product.  We  now  define, 
in  terms  of  (2.33), 


(2.36) 


A.  u  :=  l  [M  ,R  .]U  . 

S  M  J=1  M.J  m.j 


for  u,  given  by  (2.34),  in 


25  :=  (u  €  3f:  each  R  .  of  (2.34)  is  in  2)  ,  with 

y  M-J  M-'/' 


2  |M  .  R  .!?<“>}• 

s=i  m.  j  a  7 


It  is  easy  to  verify  that  this  definition  of  A^,  2)^  is  independent  of 

the  particular  choices  of  orthonormal  bases  {U  .:  j  =  1 . J (p)} 

M » J 

made  for  each  21 
M 
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Since  each  25  ,  is  dense  in  3 C  we  have  the  set  {finite  sums 

A*.'/' 

(2.34)  with  each  R  .  €  2)  ,}  dense  in  #  so  A.  is  densely  defined. 

1  '  At,  J  At.'/' 

k  ^  k  k  * 

If  u  -*  u  in  3f  (with  each  u  €  2 5^)  and  also  A^u  -»  w  in  , 

then  each  R^  -»  R  in  I  and  M  .  R^  converges  in  SC  , 

At.j  M.  j  ~M .'A  At.j  & 

Ay 

necessarily  to  M  .  R  .  .  with 
-/i.'A  p,j 

I  I.  IM  ,  R  .  |2  =  |w|2  <  « 

^J  ~At.j  M.J  * 


so  u  €  2).  with  A,u  =  w  :=  lim,  A.u 
'A  k 


I.e.,  A,  is  closed.  Next, 


suppose  u  €  2)(A^)  .  meaning  3f-continui  ty  (on  2)^)  of  the  functional  : 

u  -*  <A^u,u>  so  we  have  an  identity  <A^u,u>  =  <u,u>  for  some  u  €  %  . 

With  the  obvious  notation,  taking  u  :=  R  .  U  .  gives 

At.j  At.j 


<M  ,  R  . ,  R  =  <R  ..ft  > 

~Al.'A  ft.J  At.j  &  At.j  At.j 


for  R  .  €  2 5  , 

At.j  At.'/' 


whence,  as  M  is  self-adjoint,  we  have  R  .  €  2)  ,  and  ft 

-At.'/'  J  At.j  At.'#'  At.j 

M  /,  ,•  •  As  in  (2-35),  this  (for  each  At.j)  gives 

r* *  'r  M*  J 

I,  X.IK,  ,  3,  J2  =  lil*  <  “ 


“M  ^j  ~At.j  At.j 


whence  u  €  2)^  .  Thus  2)(A(j)  C  25Wj  =:  2)(A()  .  Since  one  obviously  has 


the  reverse  inclusion,  it  follows  that  A^  ,  is  self-adjoint.  □ 


We  remark  that  A^  ,  as  defined  in  (2.36),  has  compact  resolvent 
but  it  is  convenient  to  defer  proof  of  this  until  the  discussion  of 
spectral  analysis  of  A^  in  the  next  section.  The  final  task  of  this 
section  is  verification  of  the  consistency  of  (2.2)  (ii)  and  (iii). 


LEMMA  7  Suppose  A^  2!^  are  defined  as  in  (2.14),  directly  by 
application  of  Theorem  2,  and  also  as  in  (2.36).  Then  these  definitions 
are  equivalent 
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Since  each  2 J  ,  is  dense  in  9C  we  have  the  set  {finite  sums 
P.'P 

(2.34)  with  each  R  .  €  2)  ,}  dense  in  K  so  A,  is  densely  defined. 

M » J  M  •  t 

k  ^  k  k  * 

If  u  -»  u  in  3f  (with  each  u  €  2 )^)  and  also  A^u  -»  w  in  $  , 

k  ~  k 

then  each  R  .  -»  R  .  in  1  and  M  ,  R  .  converges  in  9 C  , 

M.J  M.J  M.j 

necessarily  to  M  ,  R  .  .  with 

(JL,  j 


so  u  €  2)  with  A.u  =  w  :=  lim,  A.u  .  I.e.,  A.  is  closed.  Next, 

\p  k 


~\p 


suppose  u  €  2!(A^)  ,  meaning  3f-continui ty  (on  2)^)  of  the  functional  : 

u  -»  <A^u,u>  so  we  have  an  identity  <A^u,u>  =  <u,u>  for  some  u  €  . 

With  the  obvious  notation,  taking  u  :=  R  .  U  .  gives 

M.j  M.J 

<M  R  R  =  <R  ..ft  4>  for  R  .  €  2)  . 

M.J  M.j  *  M.j  M.j  M-J  M-* 


whence,  as  M  is  self-adjoint,  we  have  R  .  €  25  ,  and  ft  .  = 

~M.^  M.j  M.^  M.J 

‘  AS  ln  ^2‘35^’  this  (for  each  M.j)  gives 

I„  I,IM.  ,  B.  ,l2  -  IS  I,  <  " 


-M  ~M.j  M.j 


whence  u  €  2)^  .  Thus  2D( A^)  C  2)^  =:  2J(A^)  •  Since  one  obviously  has 
the  reverse  inclusion,  it  follows  that  ,  is  self-adjoint.  □ 


We  remark  that  A^  ,  as  defined  in  (2.36),  has  compact  resolvent 
but  it  is  convenient  to  defer  proof  of  this  until  the  discussion  of 
spectral  analysis  of  A^  in  the  next  section.  The  final  task  of  this 
section  is  verification  of  the  consistency  of  (2.2)  (ii)  and  (iii). 


LEMMA  7  Suppose  A^  2)^  are  defined  as  in  (2.14),  directly  by 
application  of  Theorem  2,  and  also  as  In  (2.36).  Then  these  definitions 
are  equivalent 
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Proof-  For  u  of  the  form  (2.15)  with  U  €  we  have  u  €  25  ' (11-*  if 

r1  ' 

and  only  if  R  €  2)  ,  C  1/.  .  Since  elements  of  ^  are  smooth,  this 

p,^  0  p 

gives  u  €  II  and,  from  (2.16),  etc.,  we  have 

(A  +  *  +  X)[RU]  =  [(M  +  p  +  X)R]U  =:  ftU  +  Xu 

so  R:=  €  a  implies  RU  €  3X1(1  =  &/111)u- 
Conversely,  RU  €  means  (ft+X)U  €  3f  whence  ft  €  2£  so 

R  €  2)  , .  .  For  either  of  the  definitions  one  obtains  a  closed  operator 
fJ-.'P 

and  the  span  of  such  u  =  RU  is  dense  in  each  graph.  Thus,  the 
definitions  of  A.  ,2).  coincide  .  □ 

'W/  \h 

3  SPECTRAL  THEORY:  CONTINUITY 

We  will  be  considering  the  eigenvalues  of  A^  (defined  as  in 
Theorem  6  for  p  €  2?  with  suitable  lower  bound  ),  taken  in 
increasing  order  with  multiplicities: 

(3.1)  X^*)  $  X2(*)  $  ...  ->»  . 

as  a  sequence  of  nonlinear  functionals  of  p  .  The  principal  result  of 
this  section,  after  verifying  (3.1),  is  that  each  functional:  P  »  X^(^) 
is  continuous,  topologizing  ^  in  ^  (with  a  suitable  one-sided 
estimate:  p  >  <p  )  . 

From  the  proof  of  Theorem  6  (and  under  those  hypotheses)  we  know 

that  each  of  the  operators  ^  (p  =  0,p^,...)  is  self-adjoint  with 

compact  resolvent  and  so  we  have  eigenpairs^ *  {[a  ,  ,  y  ,  ]:k  = 

p,  K  p,  K 

1,2,...}  such  that,  for  each  p  =  0,  p^ . 


^That  is,  y  .  €  2)  ,  C  1  with  M  y  .  =  a  ,  y  ,  .  (We  are 

using  ct„  to  denote  the  eigenvalues  of  M  ^  to  avoid  confusion  with 

the  eigenvalues  {p^,...}  of  S  or  {X^,...}  of  A^  .  We  remark  at 

this  point  that  in  standard  Sturm-Liouville  theory  one  shows,  using 

2 

properties  of  the  initial  value  problem  for  [Mq  +pa/r  +  ^jy  =  0  ,  that 

these  eigenvalues  are  simple  (strict  inequalities  in  (3.2)  (ii)  with 
certain  nodal  properties  for  the  eigenfunctions.  For  p  as  rough  as 
v.ot-0  -i  t-  ic  tint  nlnar  that  this  remains  valid. 
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(3.2)  (i)  {y  ,:k  =1,...}  Is  an  orthonormal  basis  for  £  ; 

[J.  f  K. 


(ii)  V 1  *  CTu.2  1 


LEMMA  8  Each  of  the  functions  y  ,(r)U  .  (for  p  = 

k  =  1,2 _ ;  j  =  1 . J(m))  is  an  eigenfunction  of  icith 

corresponding  eigenvalue  o  ,  .  This  set  of  functions  is  an 

jj.  *  K 

orthonormal  basis  for  It  . 


Proof:  For  y.y  €  3C  and  U,0  €  we  have 

<yU,  y0>3f  =  <y.y>a  <U,0>^ 


so.  since  {U  :  p  =  0 _ ;  j  =  1 . J(m)}  is  an  orthonormal  basis 

M  •  J 

for  II  and  each  {y  ,  :  k  =  1,...}  is  an  orthonormal  basis  for  3C  ,  it 

fj. ,  K 

follows  that  {y  U  .}  is  an  orthonormal  basis  for  If  =  .  Our 

/I ,  K  fl  ,  J 

construction  of  gives 

VyU)  :=  y)U  for  U  €  ^  ,y  €  9^  : 


^yp,k  UM.  =  (V*  yM.k^UM.  J  ^M.k  y/i.k^,  j 


so  each  (y^  k  j)  is  an  eigenfunction  of  A^  .  □ 

LEMMA  9  For  any  X  €  IR  there  are  only  finitely  many  p  €  a(L)  for 

which  o( M  ,)  fl  (-“.Al  is  nonempty  so  (counting  multiplicities  in 
V/N^/lX>'p 

cr(A  ) )  the  set  { a  ,  :  a  ,  £  £}  is  finite. 

v~v/r  1  jJ,K  fl,  K 
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Proof ;  Suppose  a  £  X.  is  an  eigenvalue  of  M  .  with  corresponding 
eigenfunction  x  €  2!^  ^  C  1  ,  normalized  so  |x|^  =  1  .  Then 

a  =  <M  ,  x,x> 

2 

=  <M  x,x>  +  <<px,x>  +  <^-<p,x  > 

>  <M  x , x>  -  ll<pll  llx2llr 
-  *  [d] 

where  we  take  ll<#>ll  in  i^v)  •  From  (2.32)  and  (2.31), 

cr  >  a  Ixl2  +  (a/i  -  pQ)  -  ll<pll(e|x|2  +  CJ  . 

12 

Choosing  e  '■=  a  /ll<pll  ,  this  gives 

(3.3)  ft  i  [a  +  PQ  +  Ml  C£]/a  . 

Since  a,a,| 3^,  ll<pll,C&  are  independent  of  CT.fi,  we  see  that  a  bound  X 
on  cr  bounds  p  .  Since  we  only  consider  fi  e  ct(L)  =  (O.fi^  ,...}  , 
this  restricts  us  to  a  finite  set. 

In  particular,  given  a  ,  the  set  ^(ct)  :=  {/i  e  ct( S):  a  e  cr(M  )} 

fi  ,  T 

is  finite.  Further,  a  occurs  with  finite  multiplicity  for  each 

p  e  j«(ct)  .  For  each  occurrence  of  ct  in  [<7(14^  with 
multiplicities],  say  with  eigenfunction  y  ,  it  occurs  J(/i)  times  in 
[cr(A^)  with  multiplicities]  -  with  corresponding  eigenfunctions 
{yU^  :  j  =  l,...,J(fi)}  .  Thus,  the  set  {y^  contains  exactly 


note  that  C&  here,  coming  from  (2.31)  in  Lemma  5,  depends  only  on  a 
and  llipll  whereas  the  C£  appearing  in  the  definition  of  pg  for 
(2.32)  comes  from  (2.25),  depending  only  on  the  relation  of  a  to  the 
a  in  (1.3). 


23 


(3.4) 


I(J(M)KM:  m  €  J((cr)}  =:  ma 


eigenfunctions  associated  with  a  as  an  eigenvalue  of  . 

From  the  above  it  follows  that  {a  ,  }  is  a  discrete  set  which, 

fj. ,  JC 

when  sorted  in  increasing  order  (with  multiplicities),  we  can  relabel  as 

(Xj.Xg.--.)  with  X^  -*  00  as  in  (3.1)  with  an  associated  orthonormal 

13 

basis  for  consisting  of  eigenfunctions  of  which  we  relabel  as 

(wk  =  wicW:k  =  1*2 - >  so  =  \Wk  ’  It:  follows  that  CCT(^,) 

with  multiplicities]  is  precisely  {X^.Xg,...}  with  the  multiplicities 

given  by  (3.4) .  □ 

An  immediate  corollary  to  the  spectral  expansion  given  by  the 
eigenpairs  {[\.wk]}: 


(3.5) 


u  =  Ik  Vk  (lu|J  =  UuJ2): 


1i  ~  ^k  k 


*j,u  =  lk(\’Jk)»k  for  u  €  *here 


<  » 


v-1 


is  that  A^  has  compact  resolvent  (since  (X-X^)  0  for  X  €  cr(Ap 

{Xj,...})  .  From  (3.5)  we  easily  obtain 


(3.6) 


^  u,u>  =  ^k^kuk  for  u  €  as  in  (3-5) 


and  note,  from  this  and  (3.1),  that  <A^  u,u>  attains  its  minimum  on 

(3.7)  (u  €  2 5  :  |u|^  =  1  with  u  1  w^  for  k  <  K) 

at,  e.g.  u  =  w^  with  the  minimum  value  X^  .  It  will  be  convenient  to 
obtain  a  slightly  different  recursive  variational  characterization  of 
Xj^  .  For  p.  e  a  (L)  we  set 


13  - 

We  canonically  take  these  to  have  the  form  w  =  yU  .  (once  we  have 

M.  J  v 

fixed  the  orthonormal  bases  (U  .}  for  each  ^  )  with  y  an 

v  M.j'  PJ  J 

eigenfunction  of  M  ,  . 
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(3.8)  y  K  :=  {u  =  yU:  y  €  <3/  with  lylg,  =  1;  0  €  {U^.}; 

u  1  w^  for  k  <  K} 

with  <3/  :=  ‘S/q  for  p  =  0  and  <S+  for  ^  . 

LEMMA  10  For  the  problem 

(3.9)  minimize  <(A+'p)u,u>  subject  to :  u  €  US/*  ^ 


the  minimum  is  attained  with  the  minimum  value 


The  minimizer  u 


Is  an  eigenfunction  of  A  which  can  be  taken  to  be  . 

/v^f/  K 

Proof:  Let  X  :=  inf  (<(A-h//)  u,u>:  u  €  j^}  ■  We  already  know 

from  (3.5)  ,(3.6)  that  X-.  is  attained  at  u  =  €  U  y  so  in  (3.9) 

we  need  only  consider  u  for  which  <(A+^)u,u>  £  .  As  in  the 

argument  for  (3.3)  we  then  have,  for  u  =  yU  €  U  ^  , 

p  p,K 

\  i  <A+4»)u,u>  =  <(M^+'/')y,y>a; 

>  alyl?  +  (04  -  P0)  -  IMKelyl?  +  c6)  . 

(a/2a)  |y  +  pi  i  (^tlMCJ/a 


where  we  have  here  taken  e  =  a/2!l<pll  for  use  in  (2.31). 

If  we  consider  a  minimizing  sequence  for  (3.9),  we  see  that  only 

finitely  many  p  €  ct(S)  need  be  considered.  As  there  are  then  only 

finitely  many  relevant  (U  .}  ,we  may  extract  a  (minimising) 

Mi  J 

subsequence  of  the  form  y,U  with  U  =  U  .  (pt.j  fixed),  y,  €  3/  ,  and 

it  J  h 

<(A+'/')uk,uk>  =  <(?JM-W')yk.yk>a  -*  ^  • 

Further,  we  have  { lyk 1^}  bounded  so  we  may  also  assume  y^-*  y  (weak 
convergence  in  *3/)  . 
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Choosing  X  large  enough  that  (M^+^h-A)  :  H  -*  V  is  strictly 
monotone,  the  quadratic  form  [y  -*  ^(M^+^+-X)y ,y>g,]  is  convex  and  so 
lower  semi  continuous  with  respect  to  weak  convergence  in  *5/  .  Thus 

<(M  / W0y.y>  =  <(M  -WH-X)y,y>  -  X 

r*  r4 

i  lim  inf  <(M^+^+X)yk>yk>  -  X 
=  lim  inf  ^M^y^y^  =  X 

and  the  minimization  (3.9)  is  attained  at  u  =  yU;  set 

z  :=  (M/'/Oy  €  ■ 

Let  =  {y  €  X:  yU  €  (Wj . w^}}  with  ®  :=  ^  j 

for  the  fixed  (pi, j)  .  Then  yU  1  wk(K  <  K)  in  3f  precisely  when 
y  1  yk(k=l,  . . . ,K’ )  in  X.  .  We  set 

:=  {y  €  <y:  lylj  =  1  and  y  1  yk(k=l . K’)}  , 

p(t)  =  p(t;y)  :=  <(M+^X) (y+ty)  ,y+ty>a  for  y  €  V* 

=  (X+X)  +  2t[<z.y>a  +X  <y.y>a]  +  t2  <(M^+^X)y,y>a 

with  X  as  above.  The  minimization  property  of  y  ensures  that  p(t) 

is  minimized  at  0  for  any  y  €  9*  with  y  1  y  .  Thus,  <z,y>^  =  0 

—  ^ 
for  such  y  .  Since  z  is  orthogonal  (in  the  sense  of  the 

duality  corresponding  to  the  &  inner  product)  to  everything  in  y 

which  is  orthogonal  to  y  ,  i.e., 

z  1  (y  €  <H:  y  1  sp  (^ . yR.  ,y}>  . 

we  must  have  z  €  sp{y^ . y^,  ,y)  .  Hence  z  €  <3/  C  so  y  €  2^  ^  and 

z  =  M  ,  y  .  Also,  for  k  =  1.....K’  we  have 
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<^k,Z>3C  "  y>® 

=  K^\  =  \'<K’^  =  0 

where  >^,(k‘<K)  is  the  eigenvalue  corresponding  to  w^,  =  y^U  . 

_  —  —  nw/  —  »V  a 

Hence,  z  £  sp{y},  i.e.  z=Xy  for  some  A  .  Clearly  A  =  X  i  A^.  and 
u  : =  yU  is  an  eigenfunction  of  with  A^u  =  (M  ^  y)  U  = 

zU  =  A  u  .  The  ordering  (3.1),  i.e.  the  definition  of  A^  ,  then 

<v*  —  — 

ensures  A  >  A^  so  A^u  =  A^u  .  To  within  the  arbitrariness  in  the 
specification  of  the  eigenfunctions  we  can  take  w^  =  u  .  □ 

This  argument  is  essentially  the  Courant  Minimax  Theorem  (cf.,  e.g. 
[6]),  adapted  to  the  present  definition  of  A^  .  The  characterization 
by  (3.9)  permits  us,  as  in  [5],  to  show  the  continuous  dependence  of 
each  eigenvalue  A^  =  A^(v^)  (and  of  each  corresponding  eigenfunction 
w^  ,  to  within  the  arbitrariness  inherent  in  specif icai ton  of  the 
eigenfunctions)  on  ,  topologized  by  the  9  norm. 

_  _  ^ 

THEOREM  11  Let  \p  =  -*  \p  strongly  in  9  and.  assume  there  is  a 

bounded  sequence  (<p  =  <pi)  in  (<Sfy)  such  that  yp^  >  .  Then  as 

i  -*  00  one  has 

(3.10)  \  = 

for  each  k  =  1,2,...  .  Correspondingly ,  we  have 

(3.11)  wk  =  wkii  :=  wk(^.)  -»  wk  :=  wk(^)  in  II  =  T  (0) 

to  within  the  arbi trariness  inherently  associated  with  our  specification 
of  the  eigenfunctions . 


Proof:  The  argument  is  essentially  the  same  as  the  corresponding 
argument  in  [5],  inductively  using  the  variational  characterization: 
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(3.12)  XK(^)  =  min  (<(A+v//)u,u>:  u  €  U^  K(*)}; 

^ K(*)  :=  (U  =  yU:  y  e  <y  with  ly  =  1;  0  =  0^.; 
u  1  w^(^)  for  k  <  K} ; 

=  arg  min  (<(A+^)w.w>:  w  £ 

given  by  Lemma  10.  The  inductive  hypothesis  is  to  assume  the  result 
(3.10),  (3.11)  known  for  k  <  K  and  fix  K  .  We  now  write  X  =  X^, 
w  =  w^ ,  X,  and  w  for  X^sX^^.w^sw^^.X^,  aad  w^  . 

We  first  wish  to  show  that  lim  sup  X^  £  X  .  To  this  end,  obtain 

<v  <v  — 

u  =  u.  by  applying  the  Gram-Schmidt  procedure  to  (w1  . 1  . ,w) 

1  1 , 1  K.— 1 , 1 

so 


(3.13) 


"  *  N("  '  4<K  Sc  wk> 


where  N  =  N.  is  a  normalizing  constant  and,  noting  the  orthonormality 


of  {w^  =  wk  ^.k  <  K)  ,  we  have  i  =  N  <  w^.w  >  .  Since,  by 

inductive  hypothesis,  w^  -*  w^  and  wk  1  w  for  k  <  K  ,  we  have 


the 


(3.14) 


N  =  1^1  -»  1 ,  Ck  =  Cki'*0  for  k<K- 


From  (3.13),  (3.14)  it  follows  that  u  -»  w  in  If  .  Actually,  we  know 

that  w  has  the  form  yU  with  U  =  U  .  (some  fixed  u,  j)  and  y  €  <5/ 

pt.J  ' 

with  [yl^  =  1  and  -  y  =  X  y  .  In  (3.13)  we  have  =  0  for  any 
w^  not  corresponding  to  the  same  (m.j)  so  we  can  set 

y  :=  Ny  -  lk<K  °k  yk 


where  w,  =  y,  U  .  (any  pi.j)  and  have  u  =  y  U  (same  U  as  for  w) 
K  K  /i ,  J 


~  2, 


We  have  y  =  y .  -*  y  in  V  whence  also  (y }  is  convergent  in  fP 
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Now  Lemma  10  gives 


X  "  XK,i  -  <(£+'^u>  u  > 


=  <(MM+^)y-y>a;  +  <'^.y2>a 


We  have 


<(M  +^)y.y>  -*  <(M  +^)y,y>  =  X 

r*"  r* 


as  y  -*  y  in  “3/  and 


<  'P  -  ^  ,  y  >-»o 


as  y^-»  y  in  9s  and  'P  -*  in  2P  (even  weak  convergence  would 


suffice).  Thus,  lim  sup  X  £  X  . 

We  now  wish  to  show,  conversely,  that  lim  inf  X„  .  >  X  ,  giving 

it,  i 

(3.10),  and  that  (3.11)  holds.  Each  w  =  wv  .  has  the  form  yU  by  our 

it,  i 

specifications  and,  as  for  Lemma  10,  the  upper  bound  on  X  which  we 

have  just  obtained  restricts  attention  to  U  =  U  for  a  finite  set  of 

J 

relevant  (/i,j)  .  Thus,  possibly  subdividing  (w^  into 

14  - 

alternative  subsequences,  we  may  assume  a  fixed  U  and  that  this  U 

is  to  be  used  in  specifying  w  .  For  each  w  =  w„  .  ,  then,  we  have 

I\,  1 

w  =  yU  with  y  €  <3/  ,  ^1^  =  1  ,  and  ^  y  =  Xy  .  The  same  estimate 

as  in  Lemma  10  (recalling  the  assumed  boundedness  of  {<*>  =  W^})  bounds 
(y  =  y.}  in  *3/  so  we  may  assume  (again  possibly  taking  a  subsequence) 
that  {y.}  converges  (weakly  in  *3/  so  strongly  in  3  )  to  some  y  €  °y 

with  lylg,  =  1  .  For  k  <  K  we  need  consider  only  (y^ . y^,}  = 

(y:w,=yU  with  k  <  K)  as  earlier,  except  that  y.  =  y.  .  now  (but  we 
are  considering  the  fixed  U  as  for  w,w)  and,  similarly,  y.  :=  y.  (>//); 

J  J 

note  that  the  corresponding  indices  are  independent  of  i  by  our 
inductive  assumption.  Then 


^^These  alternatives  would  correspond  to  equally  valid  ways  of 
specifying  w^  ,  as  is  shown  by  the  subsequent  argument. 
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<*•  V* =  <y’L '  V*  - 0 


since 

=  o 

and,  inductively, 

(3.11)  giving 

w.  =  y  .  -»  w,  =  y  . 
k  k 

in 

* 

corresponds 

t°  lyj  -  yj  Ij  -»  o 

.  Hence,  in  the  limit 

<y. 

>~  =  0  and 

at 

<yU,wk>#  =0  for 

each  k  <  K  . 

We  also  have 

y2>  0 


-  is  *  2 

since  we  have  assumed  \p  -*  ^  strongly  in  &  and  {y  }  is  bounded 
in  3s  by  Lemma  5.  Now  choose  X  large  enough  that  (M^+  4>  +  X) :  3/  -* 

3/  is  (strictly)  monotone  so  the  functional:  y  »  <(M  +  4>  +  X)y,y>  is 

(strictly)  convex  on  ‘S/  and  so  lower  semicontinuous  with  respect  to  the 
weak  topology  of  3/  .  We  have 


x  =  <^~M+  *  +  x)y,y>a  -  ^  +  <xP  -  'P-  y2>a 
so,  as  y  -»  y  in  3/  ,  we  have 


lim  inf  X  £  <(M  +  'P  +  X)y,y>  ”  A 

=  <(M  +  ^)y.y>a  =  <(A  +  ^)(yu) .yu>x 

>  min  (<(A  +  ^)u,u>:  u  €  ^  . 

This  shows  that  X  =  X^(^)  ^  along  subsequences  for  which 

wv  .  =  y.U  with  U  fixed  and  (y.)  weakly  convergent  in  3/  .  The 
uniqueness 


This  corrects  a  minor  error  in  [5]  where,  at  the  corresponding  point, 
only  weak  convergence  in  &  was  assumed  -  which  seems 

inadequate  if  one  has  only  weak  convergence:  y.  -*y  in  3/  .  Note  that 

-  2 

we  need  only  lim  sup  -  ^,y  >  i  0  so  weak  convergence  would  be 

adequate  if  supplemented  by  a  one-sided  bound:  \p  <  with  strong 

converence:  0^  -»  0  in  S5  '  It  remains  open  as  to  whether  weak 
convergence  could  suffice  in  general. 
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This  completes  the  inductive  step  and,  since  the  inductive 
hypothesis  is  vacuous  for  K  =  1  ,  the  proof  of  the  Theorem  is  complete 
by  induction  on  K  .  D 

4  THE  APPROXIMATION  SCHEME 

The  method  of  generalized  Interpolation  [3],  [4],  is  a  quite 
general  approach  to  the  approximate  solution  of  ill-posed  problems. 
Typically,  one  must  first  observe  the  equivalence  of  the  problem  to 
specification  of  the  values  for  a  sequence  of  functionals  (A^*)}  ~  but 
here,  as  in  [5],  the  nature  of  the  problem  already  presents  it  in  this 
form. 

The  simplest  version  of  the  method  is  the  procedure  (P^) 

18 

described  in  the  Introduction.  The  relevant  hypotheses  are: 

(4.1)  The  norm  IMI  (determining  a  reflexive  Banach  space 

)  topologizing  the  relevant  potentials  is  such  that  if 
'p  -*  vj)  weakly  in  9  with  llvp  II  -*  ll^ll  .  then  \p  -*  \p 

V  "  l)  "  ”  u 

strongly  in  9*  . 

(4.2)  Weak  convergence  yp  yp  in  3^  implies  ^('iO 

for  each  k  =  1,2 . 

(4.3)  The  constraint  set  ♦  C  9  is  such  that  the  problem: 
A^.f'/')  =  A,  for  k  =  1 ,2,  .  . .  with  yp  € 

has  at  most  one  (minimum  norm)  solution  yp  . 

Under  these  hypotheses  (4.1)-(4.3)  it  is  a  general  result  [3]  that 


18 

The  property  (4.1)  is  referred  to  as  the  "Ef imov-Steckin  property”. 

We  refer  to  (4.2), (4.3)  briefly  as  "(weak)  continuity"  and  "uniqueness” , 
respectively. 
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(4.4)  ^  -*■  p  (strongly  in  ^  )  as  N  -*  00 

where  we  assume  the  data  {\^}  is  consistent  (i.e.  a  solution  exists  in 

(4.3))  and,  as  for  (P^),  each  ^  (N  =  1,2,...)  is  defined  as  the 

minimum  norm  element  of  P*  subject  to  matching  the  given  values  X^ 

of  X^(*)  for  k  =  1,  ....  N  . 

Rather  than  prove  the  result  in  this  form,  we  turn  instead  to 

consideration  of  a  more  general  version  which  permits  the  use  of 

(implementable)  approximate  procedures  for  the  ’N-th  stage' 

computations.  Before  doing  this  we  comment  on  the  hypotheses. 

It  is  known  (cf.  e.g.  [2])  that  (4.1)  holds  for  any  uniformly 

convex  Banach  space,  in  particular  for  Hilbert  spaces.  Our  major 

effort,  to  this  point,  has  been  to  show  that  one  obtains  continuity  of 

the  eigenvalues,  viewed  as  nonlinear  functionals  on  the  (radial) 

potential  p  ,  using  anorm  convergence  in  the  specific  space  $  and 

subject  to  a  lower  bound  condition.  Our  first  observation  is  that 

(4.1),  (4.2)  need  only  hold  on  the  constraint  set  P  .  We  will 
19 

assume  : 

(4.5)  The  constraint  set  P  is  in  and  for  each 

^-bounded  subset  ^  C  ^  ,  there  exists  a  suitable  v 
and  a  number  m  such  that  each  p  £  p^  has  a  lower 
bound  <p  €  ($  )*  with  P  2  V  .  ll<#>H  i  m  . 

and  obtain  the  (restricted  form  of  the)  condition  (4.2)  by  requiring 
compact  embedding  :  9s  -*  &  .  We  would  like  to  permit  consideration  of 

potentials  p  involving  (radial)  measures  and  note  that  our  efforts  in 
working  with  such  a  weak  space  as  $  do,  indeed,  have  the  value  of 

19 

The  simplest  form  of  this,  of  course,  would  be  to  have  p  >  0  for 
P  €  or,  slightly  more  generally,  a  one-sided  condition  that  p  >  C 

for  a  constant  C  depending  only  on  the  SP  norm  of  P  . 
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permitting  this,  even  after  the  norm  is  strengthened  (defining  2?^  )  to 

have  this  compact  embedding.  Note  that  we  are  not  assuming  that  $ 

20  * 

itself  is  compact  in  2?  but  only  a  relative  pre-compactness  in 
of  sets  bounded  with  respect  to  2^-norm  without  having  to  specify  any 
particular  such  #  bound  in  specifying  ^  . 

The  uniqueness  property  (4.3)  is,  at  present,  terra  incognita  for 
(EVP),  even  for  the  case  of  radial  potentials.  In  the  one-dimensional 
case  (fi:=  (-1,1)  C  IR^)  radial ity  just  means  that  the  potential  is 
known  to  be  symmetric  on  the  interval  and  that  is  known  to  ensure 

21 

uniqueness  [3].  This  suggests  the  possibility  that  (4.3)  may  hold  for 
quite  general  C  9  ,  satisfying  (4.5),  but  this  remains  entirely 

conjectural  at  present.  Here  we  take  the  uniqueness  condition  (4.3)  as 
an  a  priori  hypothesis  without  investigating  specific  settings  (i.e. 
more  concrete  conditions)  permitting  its  direct  verification. 


This  assumption  (corresponding,  e.g.  to  an  assumed  a  priori  bound  on 
'P^  in  a  space  as  SP  )  would  permit  a  simpler  approach.  The  map 

A  :  yp  »  OjM - ]:  -*  K°° 

00 

(taking  !R  with  the  product  topology)  would  be  a  continuous  injective 
map  from  a  compact  Hausdorff  space.  By  a  standard  result  of  point-set 
topology,  A  would  then  have  compact  range  and  a  uniformly  continuous 
inverse.  The  (uniform)  continuity  of  the  inverse  would  mean  that,  in 
specifying  'P  ,  we  have  assumed  away  the  i 1 1-posedness  of  the  inverse 
problem  (EVP).  The  difficulty  lies  in  justification  of  any  a  priori 
bound  on  the  potential  yp  . 

21 

An  interesting  stronger  conjecture  is  that  one  might  be  able  to 
recover  yp  from  knowledge  only  of  thoseeigenvalues  of  associated 

with  purely  radial  eigenfunctions,  i.e.  from  o(Mq  ^)  •  This  seems 

unlikely,  however,  as  in  the  one-dimensional  case  it  would  correspond  to 
knowing  yp  symmetric  but  only  giving  alternate  eigenvalues  -  those  with 
even  eigenfunctions. 
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We  turn  now  to  the  more  general  approximation  procedure,  relaxing 
(P^)  somewhat.  For  this  we  assume  that  we  are  given  6^  >  0  and 
positive  functions  ^(*3  ^  0  for  k  €  N  and  .  The 

procedure  (at  this  ‘N-th  stage')  is  then: 

(pa.„)  {+  £  £„.*(+)  lor 

k  =  1 . N)  and  select  ^  such  that  11*^11*  <,  inf 

(IWI*:  *  €  $N}  +  5n  • 

We  note  that  (P&  does  not  determine  ^  uniquely  . 

THEOREM  12  Let  ty*  be  as  in  Theorem.  11  and  let  9s*  be  a  reflexive 
Banach  space  ( with  norm  ll»ll  )  embedding  compactly  in  &  and 

22  it 

satisfying  (4.1).  Let  be  a  closed  convex  subset  of  ^  C  J 

satisfying  (4.5).  Assume  0  <  5^  -»  0  and  0  <  e„  .  (•)  -►  0  for  each 

fixed  k  ,  uniformly  on  -  bounded  subsets  of  .  The  operator  A ^ 

is  defined  (as  in  Section  2)  for  radial  potentials  ^  ^  and  the 

spectrum  :  k  =  1.2,...]  =  ct(A^)  is  as  in  (3.1),  in  increasing 

order  with  multiplicities.  Suppose,  for  k  =  1,2 .  we  have 

-  23  - 

A^  =  ^(*)  for  some  unique  \p  €  .  Then  (4.4)  holds  for  any 

sequence  in  **  obtained  by  the  procedure  (P&  ^)  for  N  = 

1.2,...  . 

Proof:  The  first  observation  is  that  ^  ^  so 

(4.6)  ll^ll*  <  11*11*  +  5n  ;  lim  sup  Ihf^ll*  <  11*11*  - 


22 

It  is  sufficient  that  **  be  closed  in  the  weak  topology  of  #*  . 

23  _ 

It  is  sufficient  that  *  be  unique  among  minimum  norm  solutions.  The 
general  condition  (4.3)  asserts  that  for  any  eigenvalue  sequence  (A^) 
which  can  arise  (consistency)  the  minimum  norm  solution  in  is 

unique  as  asserted. 
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We  may  (extracting  a  subsequence  if  necessary)  then  assume  ^  -*  P- 

»v 

weak  convergence  in  9*  for  some  p  .  As  we  assumed  P*  closed  and 
convex,  it  follows  that  p  €  V  and  we  will  show  that  p  is  a  (minimum 
norm)  solution  of 

(4.7)  \(p)  for  k  =  1.2.... 

<"v  - 

whence  p  =  p  by  the  assumed  uniqueness.  The  uniqeness  of  the  limit 
shows  that  the  possible  extraction  of  a  subsequence  above  was  nugatory: 
one  has  ^  -*  p  for  the  full  sequence  {\//^} . 

To  see  (4.7)  for  p  ,  note  that  boundedness  of  {^}  in  9*  gives 
^  k  :=  ^.k^N^  0  so  “*  \  by  ^Pa.N^  "  considering  only 

N  £  k  ,  of  course  for  each  k  =  1,2, ...  .  On  the  other  hand,  the 
assumed  compactness  of  the  embedding:  3^  -»  9  means  that  weak 
convergence:  p^  p  in  9*  implies  strong  convergence:  p  in 

9^  .  The  condition  (4.5),  with  boundedness  in  9*  of  {^}  also  gives 
the  "lower  bound  condition"  (^  2  <p^)  of  Theorem  11,  so  Theorem  11 

for  each  k  whence  A^t^)  = 

At  this  point  we  have  weak  convergence  ^  -*•  p  in  9*  (along  the 
subsequence).  The  convexity  of  the  norm  gives  lower  semicontinuity  with 

'V 

respect  to  the  weak  topology  so  p^  -*  p  implies 

(^•8)  \\p\\M  i  lim  inf  ll^ll*  *  ll^llM  . 

Since  p  is  a  minimum  norm  solution  of  (4.7)  by  assumption,  the 
solution  p  cannot  have  smaller  norm;  hence  Ityll^  =  Ityll^  from  (4.8) 

A/  _  _ 

and  p  is  also  a  minimum  norm  solution.  The  assumed  uniqueness  of  p 
then  implies  p  =  p  .  As  noted  above,  this  gives  weak  convergence 
'/'jl  -*•  P  in  9 ^  (along  the  full  sequence)  without  yet  using  (4.1). 


or  present  purposes  (3.10)  suffices;  (3.11)  is  relevant  only  as  part 
of  the  inductive  argument  for  Theorem  11. 


applies  to  give  \(^N)  \W 
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Now  if  we  combine  (4.6)  with  (4.8)  we  see  that  NfyH*  -»  Ityll^  . 
This,  with  the  weak  convergence,  gives  (4.4)  subject  to  the  assumption 
(4.1).  □ 

For  implementation  we  note  that  one  does  not  attempt  to  construct 
'Pjj  and  need  not  even  construct  p^  €  directly  as  in  (P^  ^)  .  If 
one  could  produce  any  element  p^  €  9*  for  which  one  would  have  an 
estimate  (for  some  ^  as  in  (P  ^),  N  =  1,2,...): 

(4-9)  S  -  V»  *  W  ■ 

with  ^(*)  0  uniformly  on  9*  -  bounded  sets  in  P*  ,  then  as  an 

immediate  corollary  of  Theorem  12,  one  also  has  p  -*  p  in  9  as 
N  -»  00  .  We  will  not,  however,  attempt  to  reduce  the  proof  of 
convergence  of  our  computational  implementation  to  Theorem  12,  but 
instead  will  use  an  essentially  similar  argument  to  prove  convergence 
directly. 

Parametrized  by  h  >  0  ,  we  will  need  a  family  of  computational 

25 

approximations  ^(h)  to  P  and  an  algorithm  which  takes  h,N,  and 
(a  representation  of)  p  €  ^K(h)  as  inputs  and  returns 

Aj,(*:h)  =  l\{r.h) . y+:h)]  . 

approximating  [A,(^/) . '^('Z')]  •  Reasonable  properties  of  such  a 

computational  procedure  would  be: 


25 


One  could  attempt  a  finite  element  discretization  of  (A +p)  from 


(2.11),  using  a  finite  element  subspace  corresponding  to  a  mesh 
parameter  h  ,  if  p  were  moderately  smooth  (or  first  approximate  p 


by  a  smoother  p  )  .  In  view  of  the  analysis  above,  one  might  more 
plausibly  use  such  finite  element  discretizations  to  (M  +  p)  for 

r* 

relevant  y  ,  assuming  a(S)  accurately  known.  This  effectively 


produces  (sparse)  n(h)  x  n(h)  symmetric  matrices  whose  ’first’  N 
eigenvalues  could  be  computed  and  taken  as  giving  A^(^;h)  .  In 

general,  such  a  procedure  would  give  p;h )  -*  A ^(P) 


as  h  -»  0  . 
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(4.10)  Given  any  \p  €  V  there  exists  ^  (h)  with 

(VI  'V 

^  6('//;h)  and  given  any  ^  €  ^(h)  there  exists 
^  €  $  with  ll'/'-'/'H*  i  S('P\ h)  there  exists  >p  €  with 
||v^-^||^  £  df^h)  where  6(*;h)  -»  0  as  h  -»  0  uniformly 
on  S6^  -bounded  sets;  we  may  also  assume  (4.5)  for 
(h)  ,  uniformly  in  h  for  small  h  ; 

(4.11)  l^k(^:h)  “  ^('JO  1  ■*  0  as  h  -*  0  for  each  fixed  k  = 
1,2,...,  uniformly  on  #  -bounded  sets. 

Without  further  concern  for  the  details  of  possible  construction  of  such 
algorithms,  we  indicate  how  the  availability  of  a  computational 
implementation  satisfying  (4.10),  (4.11)  could  be  used  to  obtain  a 
computable  sequence  {^}  converging  in  3^  -norm  to  4>  . 

We  wish  to  replace  the  approximation  procedure  (P  ^)  by  a  more 
explicitly  implementable  computational  procedure: 

(Pc  jj)  Choose  h  =  h^  small  enough  that 

#N(h)  :=  &  €  **(h>  :  lV^:h)  "  V  *  ^N.k  for  k  *  N} 

is  nonempty  and  select  ^  €  ?  :=  ^(h^)  such  that 

ll'^N,l*t  *  UN  +  where  :=  inf  t11^11*  1  'P  £  • 

The  actual  computation  involved  in  (Pc  N)  would  be  the  use  of  some 
(standard)  algorithm  for  nonlinear  constrained  optimization  to  minimize 
Ityll^  (using  a  stopping  criterion  giving  approximate  minimization  to 
within  6^  >  0  of  the  infimum  u^)  subject  to  the  constraint  : 

<v  <v 

^  ^  .  The  computational  difficulty  of  this  will  depend  on  the  nature 

of  II  *11^  ,  on  the  sizes  of  ^  on  the  size  of  h^  and  the 
computational  difficulty  in  implementing  A^(*;h^),  etc. 
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THEOREM  13  Let  A  ,&*  .9*  ’{\)  •  be  as  f°r  Theorem  12  and  assume 

implementable  computational  approximations  ^(h) ,  an(  •;h)  are  available 
(for  small  h  >  0)  satisfying  (4.10),  (4.11)  .  Assume  0  <  5^  -»  0  and 

Zy  k  -»  0  for  k  =  1,2,  .  . .  as  N-*00.  Then ,  for  each  N  =  1.2, , 

one  can  choose  h  =  h^.  so  :=  ^(hjj)  ts  nonempty  (further  requiring 
that  h^  -»  0)  and  select  ^  as  in.  (Pc  .  For  any  such 

computed  sequence  {^}  we  have 

(4.12)  -*  4>  in  9^  -norm  as  N  -»  00 

where  ^  Is  given  by  (4.3)  . 

Proof:  Since  \p  €  $  one  has,  by  (4.10),  existence  of  '/'(h)  €  $  (h) 

with  ll'p  -  \//(h)  11^  £  5(\/); h)  -*  0  .  This  makes  (\£(h) :  0  <  h  <  h^} 

bounded  so  (4.11)  gives 

l\Wh);h)  -  \(4>(b))\  $  eN  k/2  for  k  =  1 . N 

for  small  enough  h  .  On  the  other  hand,  ^(h)  -»  \p  in  9  (a 
fortiori ,  in  $  )  and  we  have  assumed  (4.5)  for  ^(h)  so  Theorem  11 
applies  to  give 


l\0Kh))  -  \\  <  ^  k/2  for  k  =  1 . N  . 

Combining  these  gives  ^( h )  €  ^(h)  for  h  small  enough  so  then 

^(h)?50  .  Also  requiring  h^  <  h^  for  any  given  sequence:  0  <  h^  -*  0 

26  ~  ^  ^ 
lets  us  fix  h^  .  This  fixes  ^  :=  ^(h^)  and  we  can  find  ^  , 

approximately  minimizing  the  norm,  as  in  (P  „)  . 

C « i i 

To  show  (4.12),  we  proceed  as  in  the  proof  of  Theorem  12  .  Since 
'//(h^)  .  as  above,  is  in  for  each  N  and  <l»(h^)  4  +  ■  we  see  that 


26 


Computationally,  one  might  start  with  a  trial  hn  =  h^  and  then,  say, 
successively  halve  h^  until  one  can  obtain/compute  some  'p  €  ^(h^)  . 
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so  {ujj}  is  bounded  and,  as  in  (4.6),  we  have 
(4.13)  lim  sup  H'/'jjil*  i 

so  (extracting  a  subsequence  if  necessary)  we  may  assume  weak 
convergence:  >//  -»  \p  in  .  By  (4.11)  and  the  definitiion  of  ^ 

(noting  that  ^  -»  0)  we  have  X^O^)  ^  while  Theorem  11  gives 

It  follows  that  is  a  solution  of  (4.7)  with 
H'/'ll*  i  lim  inf  Ih/Lll  £  ll'jDll^  :  the  uniqueness  property  (4.3)  then  gives 
\p  =  \p  and  weak  convergence  \p^  -*  \p  for  the  full  sequence  .  From 
(4.13)  we  also  have  -»  llv/zll^  so,  by  (4.1),  we  have  (4.12)  as 

asserted.  □ 


\{  v  ■*  • 
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